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Abstract 



The relation between the restricted path integral approach to quantum mea- 
surement theory and the commonly accepted von Neumann wavefunction 
collapse postulate is presented. It is argued that in the limit of impulsive 
measurements the two approaches lead to the same predictions. The exam- 
ple of repeated impulsive quantum measurements of position performed on 
a harmonic oscillator is discussed in detail and the quantum nondemolition 
strategies are recovered in both the approaches. 
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I. INTRODUCTION 



Despite the impressive successes of quantum mechanics in explaining most experimental 
results about microscopic phenomena, a unique consistent quantum theory of measurement 
is still missing. In this field, an important achievement was made by von Neumann which 
postulated two ways for the evolution of the state vector: continuously, according to the 
linear Schrodinger equation, when time passes without a measurement on the system being 
performed; and discontinuously, according to probability laws, if a measurement is carried 
out (the so-called wavefunction collapse [|I|). von Neumann first defined this feature "a pe- 
culiar dual nature of the quantum mechanical procedure, which could not be satisfactorily 
explained" Q. 

Mensky proposed a phenomenological theory, based on Feynman path integral formalism 
0, that expresses the a posteriori dynamical evolution of a system, undergoing a contin- 
uous measurement, in terms of the instrumental uncertainty and the output of the meter, 
supposed to be known before doing the calculations. The effect of the measurement is in- 
troduced in the space of the paths by means of an influence functional, which restricts the 
integration to those paths that lie around the measurement result. The restricted path in- 
tegral approach has been applied to describe continuous [^] and impulsive [^] measurements 
in both linear and non-linear systems, to explain quantum Zeno effect and to verify the 
possibility of testing quantum mechanics through temporal Bell-like inequalities in bistable 
potentials f^. 

The most common objection raised against this approach is that it seems to treat the Feyn- 
man paths, in some sense, as real trajectories followed by the system. Mathematically, this 
means that the time evolution during the measurement is non-unitary, i.e. the wavefunc- 
tion looses its normalization. In fact, it has been shown that the obtained propagator 
is the same as the one associated to an effective Hamiltonian having a purely imaginary 
measurement term, which of course destroys the unitarity. 

The aim of the present paper is to show how this effect is equivalent, in the impulsive 
limit, to von Neumann collapse. Section II is devoted to derive, from the restricted path 
integral propagator for a measurement of infinitesimal duration performed on a generalized 
coordinate, an analytic formula with which von Neumann formalism can be recovered. In 
Section III the calculation for a series of impulsive measurements in a generic potential with 
discrete energy levels, only requiring the knowledge of the energy eigenstates, is developed 
and is applied - in Section IV - to the case of repeated impulsive measurements of position 
in a harmonic oscillator Finally the predictions of Mensky and von Neumann theories 
are compared and conclusions are drawn in Section V. 



^ "eine eigenartige Duplizitat des Vorgehens, die nicht geniigend erklart werden konnte" |jl] (p. 222 
in the original version; p. 417 in the English translation). 
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II. QUANTUM MODEL FOR IMPULSIVE MEASUREMENTS OF POSITION 



The path integral formulation of quantum mechanics provides a natural framework 
for handling continuous quantum measurements of position by restricting the integration, 
in the space of the trajectories x{t) in the coordinate x - assumed to be continuously 
monitored between times and r to those paths which turn out to be compatible with 
the experimental outcome a{t) within the instrumental error Aa This can be done by 
means of a weight functional w^a] [x] depressing the contribution of paths whose distance from 
the a{t) actually obtained (which is in general a continuous but not necessarily differentiable 
function of time) exceeds Aa. The propagator for the system is then written as a weighted 
integral: 

/■x{r)=x" ( i 1 

K[a]{x",T;x',0) = d[x]exp\- C{x{t),x{t),t)dt\ W[a][x]. (1) 



x(0)=x' [h Jo 

The probability distribution for the measurement output is a functional of a{t); if ij{x, 0) is 
the wavefunction representing the initial state of the system, it has the form 



I ma]irMa]irWd[ay 
where 

^|J[a]ix", r)= j K^a]{x\ T- x', 0)^(x', ^)dx'. (3) 



Its dispersion estimates the actual experimental accuracy with which it is possible to extract 
the information on the position x, also called effective uncertainty: 

Aa^g = 2 y ^ 11^°"^^^ ~ ~a{t)fdt Pia\d\a\. (4) 

where d(f) is the path which maximizes P^a]- Of course, in general Aagg > Aa. 

The most natural way to represent in this framework an impulsive measurement at the 
initial instant is to consider it as the limit for infinitesimal time intervals of a continuous 
one with constant result a(t) = a. Of course, in this case the probability distribution for 
the measurement results is a function of a. 

Alternatively, one can take the limit directly in the path integral expression (|l|). A simple 
form for the weight functional is the Gaussian one [Q: 

= exp |— K y \x{f) — a{tyydt\^ , (5) 

where the measurement coupling k in general should be taken constant, in order to ensure 
that the dynamics can be described by a semigroup 0. In the limit r 0, it turns out to 
be useful the position 
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in which Aa assumes the proper meaning of a width in the space of paths. In this way, 
regardless of the form of the Lagrangian involved, one obtains an analytical expression |]^: 

Ka{x",x') = \miK[a(t)=a]{x",T;x',0) = 

= lim/^^ , rf[x]exp / h-^- ^ \',J \dt\ = 



r~^0Jx(Q)=x' [Jo \h 2Aa'^T 



hm / d\x] exp <^ \, ^, ) = (7) 

e-^lx;^ir(x",0;x',0) = 



= e 2Aa^ — a;'). 

The third line follows from neglecting, in the limit, the Lagrangian term with respect to the 
measurement term, going as (this appears reasonable also from the physical point of 
view, because an impulsive measurement is assumed to induce a significant change in the 
state of the system during a negligible amount of time), and by applying the theorem of the 
mean value. K is the propagator in the absence of measurement. 

Let ilj{x,t) be the wavefunction which describes a system undergoing an impulsive measure- 
ment at the instant t, with result a. From Eqn. follows 

1pa{x, t+) = Wa{x)lp{x, t~), (8) 

where 

def (x— a)^ 

Wa{x) = e 2Aa^ . (9) 

Thus II ipait'^) II is the projection {i.e. the scalar product modulus) of ipi^x, t^) on the weight 
function Wa{x). So the quantity 

Pia) = ,.Jif"^'"^ = 7^ II Mt^) r (10) 

JJe-^^\^{x,t-)\^dxda v''^" 

represents - in analogy with Eqn. @ - the probability that the system is found in the state 
described by Wa{x). 

It should be noted that in the limit r — >^ 0, if Aa remains finite, k diverges as r~^: this means 
that von Neumann collapse is recovered by considering not only infinitesimal measurement 
durations, but also infinite coupling between the instrument and the system. Furthermore, 
in the limit of an extremely precise measurement (Aa 0), k diverges more rapidly than 
r~^, and one easily obtains 

hm P(a) = |^(a,r)|2 (11) 

Aa-*0 

in agreement with the conventional interpretation of the wavefunction. 
It is worth noting also that || ipa{t^) ||^ can be less than 1: this follows from the non- 
unitarity of the temporal evolution induced in Eqn. (|l|) by the measurement term. Thus, in 
the impulsive limit, that apparently arbitrary modification of the dynamics of the system 
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under monitoring is equivalent to the commonly accepted von Neumann postulate [0] of 
the discontinuity introduced by the measurement in the causal linear evolution of the state. 
For instance, the case of a perfect measurement which yields with certainty the information 
whether a particle is found within an interval of width Aa around the position a, is recovered 
by choosing for the measurement operator the form, discontinuous and therefore less realistic 
than the (|), 

wl^- oc 6{x - [a - Aa])e{\a + Aa] - x). (12) 

The restricted path integral approach to continuous measurements with a generic weight 
functional W[a] [x\ appears therefore simply as a generalization of the idea of perfect instanta- 
neous filtering in von Neumann theory of measurement being a smoothed version of it with 
finite accuracy and duration. The comparison between the results of the two approaches 
will be now performed in the case of stroboscopic measurements of position on a harmonic 
oscillator 



III. STROBOSCOPIC MEASUREMENTS OF POSITION 

A stroboscopic sequence of measurements is obtained when an observable is monitored in 
an impulsive way at some definite instants, equally spaced by a quiescent time AT in which 
no measurement is performed. Such a topic has been studied in detail |^ for characterizing 
Quantum Non Demolition p[ strategies for the measurement of the position of a quantum 
system. Particular advantages, in this field, can arise from applying the method developed 
in the preceding section and exploiting an energy eigenstates expansion, because of the 
uniformity to handle each form of the weight function Wa{x)^ for instance the naive one 
(0) - which expresses the usual representation of the measurement -, allowing analytical 
calculations for every potential [0. For a generic system having discrete energetic levels {i.e. 
H\l) = Ei\l)), an initial state can be developed in energy eigenstates: 

oo 

mto)) = T.ci\i). (13) 

1=1 

If the coordinate of the system is measured with results a„, n = 0, 1, . . . , iV at each of the 
instants tn = nAT, we get 

I^W}„.o......(4)) = (j[e-r.^''^Wa^_}j m^)), (14) 

where Wa is the multiplication operator corresponding to the weight function (^. The 
normalization constants relative to each measurement are factorized and therefore can be 
neglected in the calculations: because they simplify in the definition (^ of P{a), Aacs will 
not depend on them. 

By considering Eqn. ([13|) and by inserting a completeness 1 = I]m=i |"^)("^|, Eqn. ( |14|) 
is rewritten 
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|V'{a„}„=o,....iv(^iv)) = Ci\m){m\Wa^e '^^^^Wa^_, ■ ■ ■ Wa,e ^^"^^Wao\l) 



l,m=l 
CO / oo 



E E^^/C, |m), (15) 

m=l \l=l / 



where 



TV 

B^,{AT,Aa,{a^}) {m\w^, [Ue-^^'^^w^.-A 10- (16) 
By inserting times the identity operator 1 = I]^=i \ni){ni\ {i = 1,2, . . . , N), one gets 

oo /N-1 \ ( -Arp N >| 

= E n W^nTnZ. exp --- Y.EnA, (17) 



ni,n2,...,njv=l \i=l / *=1 



with 



dof 



+00 



W,f{Aa) = / M*(x)wa„(a;)Mj(x)da;, (18) 



in which the Ui{x) = {x\i) are the energy eigenfunctions in {Afj-representation. 

Let Pao,ai,...,aM-i{(^) ^e the probabihty that the A^*'^ measurement will give result a, when 
the results of all the previous measurements are known. Its dispersion can be now evaluated, 
in analogy with Eqn. (^, from Eqn. (|T3]) through Eqs. (|T^-(^H|). It is no more necessary to 
take the mean value on the measurement time as in (^, since t 0. The final expression 
is R: 



AalMan}n<N-i,N) =2 J 




m=l 



E -^mi("0, • • • , flAf-l, a)ci 



2 

da, (19) 



where Sat is the most probable result of the N measurement, and the explicit dependences 
of Aoeff and of B^i from AT and from Aa have been omitted. 

Summarizing, once eigenvalues and eigenstates of a generic quantum system with discrete 
energy spectrum are known, from the decomposition in eigenstates of the initial wavefunction 
it is possible in principle to infer directly the value of the effective uncertainty after an 
arbitrary sequence of stroboscopic measurements. The actual calculation is difficult due 
to the presence of multiple sums on the numerable ensemble of the energy eigenstates. A 
practical evaluation requires to approximate the sums and the integrals in Eqn. (|19|), by 
truncating them to finite values, say respectively A^max and i^MAx- The accuracy of such 
an approximation has to be checked for comparison with results already known by other 
methods, both analytical (if any) or numerical. 



IV. OPTIMAL MEASUREMENTS OF POSITION FOR A HARMONIC 

OSCILLATOR 

The harmonic potential shows its great importance in many fundamental problems in 
physics. In particular, it has been studied in the limit of infinite coupling for impulsive 
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measurements, to search for Quantum Non Demolition (QND) observables P]. These last 
are relevant for the detection of small displacements in the quantum limit of sensitivity for 
mechanical resonators used as gravitational wave antennas. In the present framework this 
problem can be formulated by noting that in general the effective uncertainty Aaes is greater 
than the instrumental error Aa, expressing the spreading of the paths due to the effect of the 
back-action of the meter on the measured system. However for certain observables - namely 
the QND ones the ratio of the two uncertainties can be reduced to unity by applying 
optimal measurement strategies, without violating the Heisenberg principle P|. In this 
section the method described before will be applied to the characterization of such QND 
strategies for a harmonic oscillator. The results obtained through the two measurement 
theories under consideration, expressed respectively by the weight functional (|D and (|), 
will finally be compared. 

For a harmonic oscillator, described by a Lagrangian of the form C = |mi;^ — ^mu'^x^, 
it is possible to characterize also a priori the QND strategies for stroboscopic measurements 
of position, by means of the commutation relation of the observable xh - in the Heisenberg 
picture - at different times 0: 

ih 

[xuit),xnit + AT)] = —sm{iuAT). (20) 

From Eqn. (^) follows that a series of impulsive measurements of position with infinite 
precision, performed every T/2 - where T = 27r/uj is the oscillation period -, will give 
perfectly predictable outcomes. Therefore one expects that, for impulsive measurements 
spaced in such a way, the inequahty Aoefr > Aa could be saturated. This can be seen also 
by looking to the periodical wavepacket reformation via causal dynamical evolution. In the 
absence of measurements, the wavefunction ip{x,to) = {xlipito)) of the state (|13|) will be, 
after half an oscillation period - in the same notations as above, specialized to the harmonic 
potential -, 

oo 

e-^f^i^-^o)J2c,ui{x) = 
1=1 

oo 

Y.e-^(^+'^>^^ciUi{x)= (21) 
1=1 

oo 

e-^Y^e-''-Qi-iyui{-x) = 
1=1 

-iij{-x,to), 

the third line following from the definition of T and from the symmetry properties of the 
energy eigenfunctions. 

Thus, every half-period, the wavepacket is reformed symmetrically with respect to the equi- 
librium position, except for an irrelevant phase factor. This means that, if ?/^(x,to) repre- 
sented the state after an impulsive, infinite-precision (Aa — 0) measurement with outcome 
ao, another measurement att = tQ+T/2 will yield with certainty the result — ao; at t = to + T 
this will be again ao. Therefore, in the case ao = ||^, the optimal QND strategy - in which 
Aaeff — >■ Aa, also for Aa > - is obtained by choosing the quiescent time AT = T/2. On 
the other hand, if ao 7^ 0, the optimality is reached either with AT = T/2 and alternated 



T 

1p {X,t = tQ + — 
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results Qn = (— l)"ao = 0, 1, 2, . . .), or with constant results a„ = Oq but AT = T. 
Indeed, this behavior has been already obtained by simulations based on the restricted 
path integral; here it will be recovered via von Neumann collapse, by means of the calcu- 
lation described in the previous section. A preliminar comparison of it with the results of 
Ref. should allow to test the degree of accuracy of the adopted approximations. To do this 
the conditions chosen there have been exactly replied: a Gaussian initial state with width 
0" = 5 (in the unit system in which h = 2m = 1) centered at the origin, an instrumental 
error Aa = 1 and a measurement time r ~ 10~^T. The investigation has been restricted to 
sequences of measurements with constant result a„ = Oq for n < N — 1. In any symmetrical 
single-well potential, a symmetrical state localized in the middle of the well will change its 
width during the dynamical evolution but will remain centered in the same position at all 
times. Therefore, starting from such an initial state and choosing oq = 0, the most probable 
results of a measurement is still a„ = for every n and any AT. With these values of the 
parameters, the calculation has been repeated, varying the value of AT, according to three 
methods: 



A) If an impulsive measurement is simulated 
duration r ^ Tc == 



I by means of a continuous one of short 
the computed quantities do not depend 
upon the duration of the measurement. With the choices made above, during the 
measurement the state remains a Gaussian one, having width 



a{t) = - 



1 



a + 



1 

2^ 



— a 



(22) 



with 



a 



imiOr cos LOrT 



/5 



imiOr 



2 2 



ih 



TmAa? 



(23) 



After the end of the measurement, the state evolves causally, still preserving the Gaus- 
sian form with 



a{t + T) 



a' + 
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4(2^ 



a' 



(24) 



the coefficients refer now to the unmeasured dynamics: 

imuj cos ujt „, imuj 



a 



2^ sin cut 



sin ujt. 



(25) 



By alternating these two types of evolution and by iterating this procedure N times, 
one can obtain the width of the state after stroboscopic measurements. In the 
impulsive regime, Aocfr is then easily calculated, according to its physical meaning 0: 



lira AacfT(7-) 



Aa2 



cr 



N- 



(26) 
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B) As shown in Ref. [Q, the evaluation of the path integral (|T]) can be overcome by 
writing an effective Schrodinger equation which takes into account the effect of the 
measurement through an imaginary potential term: 

Eqn. (^) can be resolved with usual numerical techniques on a space-time lattice. 
By turning on and off the measurement potential with periodicity AT, one can compute 
the values of ip^x) on the chosen lattice after the desired measuring sequence; Aoeff is 
then calculated via Eqs. (|D and (|). 

C) Method A applies only to Gaussian weight functionals. On the other hand, Eqn. ( |T9|) 
do not need this requirement; furthermore, it gives directly the value of Aqcs after an 
arbitrary series of measurements, provided that convenient approximations have been 
made. The results shown in this paper have been obtained by restricting the extremes 
of integration in Eqn. (|19D to ±10Aaes, and by truncating the sums to A^max = 20. 
After comparison with the other approaches, this method will be applied to forms of 
the weight function different from the Gaussian one 

Fig. 1 shows the typical evolution of the dispersion Aoefr of the probability P{a), as a con- 
sequence of repeated measurements with different periodicity (the cases of AT/T = \, ^ 
and I are shown). Suddenly after the first measurement, performed at the initial instant 
of time, Attes approximates the width of the initial state; to be more precise, its value is 
given by Eqn. (^): in the considered case it is Aa^s — 5.099, as evidenced by the point 
n = 1 of the graph. Due to the following measurements, the effective uncertainty reaches 
an asymptotic value Aa°g that depends on the quiescent time AT. The values of Aa^| 
given by the different methods of calculation show a good agreement: so the more flexible 
approach C can be held for tested and can finally be applied to the comparison between the 
two measurement theories of Mensky and of von Neumann, which is the aim of the present 
paper. This is shown in the last figures. 

Fig. 2 depicts the dependence of the asymptotic probability distribution for the measurement 
outcomes, Pas(o), on the quiescent time. In other words, it is considered the probability 
distribution for the result a of the A^**^ measurement (A^ = 16 has been chosen, i.e. in the 
asymptotic region as shown in Fig. 1): in the graph are visualized together many curves 
relative to different sequences of measurements, each one characterized by a quiescent time 
AT, as written on the y-axis. It can be seen with great evidence how the cases with AT 
integer or semi-integer multiple of the oscillation period are characterized by minimal dis- 
persion in the Pas{o). 

These results are summarized in Fig. 3, which shows the behavior of the asymptotic Aa^^ 
versus AT. It is obtained a periodic behavior, with minima each |^ in which Aa^l approxi- 
mates the value Aa, which identify the optimal QND strategies. The same pattern would be 
obtained by choosing oq 7^ 0, although in this case, in general, one has also a„ 7^ unless 
AT/T = i G N. As explicitly shown above, the latter condition follows from the fact that, 
starting from a state not centered in the origin, one must wait at least an oscillation period 
for the wavepacket to be reformed at the same position. Thus the definition of Aceff is no 
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longer the same for all AT, except for AT = T, 2T, . . . : for these quiescent times the QND 
strategies are recovered. 

Finally, the basical aim of Fig. 3 is the comparison between the results obtained when the 
measurement is simulated respectively by a Gaussian weight (as for Fig. 2) or by a double- 
step perfect filter as w^'^' (0). The most significant difference between the two profiles 
consists in the presence, in the case d la von Neumann, of oscillations and flexes, probably 
attributable to the discontinuity of w^'^'(x). But the essential qualitative features {i.e. the 
ones that can find experimental application does not depend upon the detailed choice 
of the weight Wa- Thus, for impulsive measurements, von Neumann wavefunction collapse 
is simply a particular case of the restricted path integral theory. 



V. CONCLUSIONS 

The restricted path integral approach to quantum measurement can be straightforwardly 
specialized to impulsive measurements, providing analytical formulae for the evolution of the 
measured system. The model leads to a filtered state centered around the measurement re- 
sult through instantaneous collapse of the wavefunction, having as width the instrumental 
error. Von Neumann ideal collapse is recovered by choosing a double-step ideal filter as 
measurement weight functional. The probability distribution of the possible measurement 
outputs has a dispersion whose behavior does not depend upon the form of the filter. In 
particular, for stroboscopic measurements of position on a harmonic oscillator, the optimal 
measuring strategies can be characterized both using an ideal or a Gaussian filter. On the 
other hand a continuous, smooth weight appears more adequate to represent in simple terms 
a physical process as the measurement. Indeed, the instrument being constituted by a col- 
lection of microscopic objects immersed in a thermodynamical environment |jlO[, it is more 



likely to give rise to a Gaussian probability distribution for the measured physical quantities. 
A Gaussian weight is also favorable from the analytical point of view, and guarantees some 
natural group properties of the evolution under continuous monitoring ||^. 
Summarizing the obtained results, the restricted path integral theory of quantum measure- 
ment reveals to be a simple and natural extension of the naive von Neumann instantaneous 
collapse theory to the more realistic case of continuous and non-ideal measurements, and 
appears as a more subtle tool for handling new fundamental problems (see for instance 
Refs. M). 



ACKNOWLEDGMENTS 

I am grateful to R. Onofrio for encouraging me to write this paper, and to him and 
C. Presilla for a critical reading of the manuscript. 



10 



REFERENCES 



[1] J. V. von Neumann, Mathematische Grundlagen der Quantenmechanik, (Springer, 
Berlin, 1932); English translation by R. T. Beyer, Mathematical Foundations of Quan- 
tum Mechanics (Princeton University Press, Princeton, 1955). 

[2] M. B. Mensky, Continuous Quantum Measurements and Path Integrals (lOP Publishers, 
Bristol and Philadelphia, 1993). 

[3] R. P. Feynman and H. R. Hibbs, Quantum Mechanics and Path Integrals (McGraw-Hill, 
New York, 1955). 

[4] M. B. Mensky, R. Onofrio, and C. Presilla, Phys. Lett. A161, 236 (1991). 

[5] M. B. Mensky, R. Onofrio, and C. Presilla, Phys. Rev. Lett. 70, 2828 (1993). 

[6] R. Onofrio, C. Presilla, and U. Tambini, Phys. Lett. A183, 135 (1993); U. Tambini, C. 
Presilla, and R. Onofrio, Phys. Rev. A51, 967 (1995). 

[7] T. Calarco, R. Onofrio, Macrorealism, Non-Invasivity and Quantum Mechanics: a quan- 
titative approach, in Proc. Conference on Foundations of Quantum Mechanics, Lecce 
1993, (Kluwer, Dordrecht, 1994); Phys. Lett. A198, 279 (1995). 

[8] C. V. Caves, K. S. Thorne, R. W. Drever, V. Sandberg, and M. Zimmermann, Rev. 
Mod. Phys. 52, 341 (1980). 

[9] W. H. Press, Numerical Recipes: the Art ofScientiGc Computing (Cambridge University 
Press, Cambridge, 1986). 
[10] A. Konetchnyi, M. B. Mensky, V. Namiot, Phys. Lett. A177, 283 (1993). 



11 



FIGURES 



FIG. 1. Evolution of Aogfj versus the number of measurements n, for some values of the quies- 
cent time AT. The continuous line refers to the analytical calculation done in Ref. |p (crosses are 
numerical results) by approximating each impulsive measurement with a continuous one of short 
duration r ~ lO^^T. The dotted line is deduced from Eqn. ([T9|), i.e. in the limit for r — > 0. The 
agreement, besides slight deviations in the initial transient, is within 1% in the asymptotic region. 

FIG. 2. Asymptotic probability distribution Pas{o) for the measurement outcomes a as function 
of the quiescent time AT/T. The curve reaches its maxima (and correspondingly the minimum 
dispersion in a) at AT multiples of half an oscillation period. On top of the graph is superimposed 
a contour plot which shows the extremal regions. 

FIG. 3. Normalized asymptotic effective uncertainty Aa^g/Aa versus the quiescent time AT/T 
for stroboscopic measurements performed with Gaussian (|5|) or perfect (^) filtering. The curves 
does not differ qualitatively and coincide at the minima, as theoretically expected ||5|. 
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